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Abstract 

The growth of matter perturbations in the f{R) model proposed by Starobinsky is studied 
in this paper. Three different parametric forms of the growth index are considered respectively 
and constraints on the model are obtained at both the la and 2a confidence levels, by using the 
current observational data for the growth factor. It is found, for all the three parametric forms 
of the growth index examined, that the Starobinsky model is consistent with the observations 
only at the 2a confidence level. 
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I. INTRODUCTION 



The present cosmic accelerating expansion jl-|29| is one of the key challenges in fun- 
damental physics and cosmology. There are basically two kinds of options to explain this 



mysterious acceleration. One is the well known dark energy an energy compo- 

nent, which has a sufficient negative pressure to induce a late-time accelerated expansion; 
the other is the modified gravity, which originates from the idea that our understanding 
of gravity is incorrect in the cosmic scale and general relativity needs to be modified. One 
of the popular modified gravities is the f{R) theory (see 36l439l| for a review), where R 
is the Ricci scalar and f{R) is an arbitrary function of R. For an f{R) model, its action 
takes the form 



S 



d^x^ —g 



-f{R) + Lr, 



(1) 



where g is the trace of the metric g^^, Gat is a bare Newton gravity constant and is the 
Lagrangian of matter. Considering a spatially fiat Friedman-Lemaitre-Robertson- Walker 
universe, whose metric is ds'^ = —dt^ + a?'{t)d'x^, and varying the above action with respect 
to gfj,y, one can obtain 



-2FH 



StiGn {Pm + Prad) + ^(Fi? - /) - 3HF 



SttGn ( Pm + -prad ] + F-HF 



(2) 
(3) 



where R = 6{2H^ + H), an over-dot stands for a derivative with respect to the cosmic 
time t, H = ^ is the Hubble parameter and F — '^-^ 



Originally, Capozziello 



dR 



40| proposed an f{R) model, f{R) = R-a/R"^ (a > 0, m > 0), 



to explain the present accelerating expansion. However, this mode 



problems, which are related to the solar-system constraints 4l|, the instabilities 42 1 



was plagued with some 



a viable cosmic evolution history with an accelerating expansion 43|] and a standard 
matter-dominated stage {4^. The main reason this model does not work is that f^Rn = 
d'^f/dR^ < 0, which gives a negative mass squared for the scalaron field. Soon, the 
aforementioned problems were solved, for example, the instabilities and the inconsistence 



with the solar-system constraint were solved in Refs. 



45 



46j, and t 



dominance was solved in Refs. 47H49| . Later, Amendola et al. 54] gave the conditions 



le problem of matter 



to obtain a viable f{R) model. Some models satisfying these conditions, the Starobinsky 



model, for an example, 
there are some models 



lave been proposed 



45 



50 



5J-|62| . Moreover, it is interesting to note that 



53( 1 in f{R) gravity, which can not only explain the present 



accelerating expansion successfully, but, at the same time, can also yield an inflation in 
the early era of our universe without a scalar field. 

Let us note that both the dark energy and f{R) gravity can explain the present ac- 
celerating expansion. However, although different models can give the same late time 
expansion, they may produce different growths of mat ter p erturbations 63|. Thus, the 
studies of the linear growth of matter perturbations 64l4lOO[ | provide a particular method 
to discriminate different models. Defining the growth function 6{z) = 6p rn./ Prp (pm is the 



din (5 



energy density of matter) and the growth factor / _ ^y^^, 
that / can be parameterized as 



the authors in 



101 



102| found 



(4) 



where 7 is called the growth index and Qm is the fractional energy density of matter. If 7 
is treated as a constant, its theoretical value can be obtained by expanding the equation 
of 7 around Qm — !> which is a good approximation at the high redshift. Then different 



models lead to different theoretical va. 
for ACDM model and 70, 



11/16 [sd 



80 



82| 



ues of 7 [80|-|99|, for example, 700 — 6/11 
8l| for flat DGP model. Therefore, it is possible 
to distinguish them. By comparing the theoretical value of 7 with the observed one, one 
can hopefully single out the model which is consistent with the observations. 

However, the growth index is, in general, a function of redshift. Some works have 
been done on the evolutionary form of 7(2;). In Refs. 89N96| . the authors studied 7(z) 
with a linear expansion, 7 ^ 70 + 7qZ, and found that this form gives a very good 
approximation at the low redshift z < 0.5 and for different models 7o is different. Thus, an 
accurate measurement of 7g could provide another characteristic discriminative signature 



to discriminate different models. In Refs 



100|, we proposed a parametrization 7(2;) 



7o + 71-2/ (1 + -2), and obtained that, for wCDM and DGP models, this form approximates 
the growth factor / very well both at the low and high redshift regions. 

In this paper, we aim to examine the growth factor of matter perturbations in f{R) 
gravity and we take the Starobinsky f{R) model as an example. Let us note that the 
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density perturbations of the Starobinsky f{R) model have been studied systematically in 
the literature [91-95|. But what we plan to do here is to examine different parametric 
forms of growth index and study the observational constraints from the growth factor 
data. 



II. THE STAROBINSKY'S MODEL 



The Starobinsky's model has the form: 

f{R) = R + XsRo 



- 1 



(5) 



where and n > are two positive constants, and Rq corresponds essentially to the 



91- 



present value of the Ricci scalar R. This model has been studied in the literature [62 
95| and it has been found that, when n > 2, all known the laboratory and Solar system 
tests of gravity can be satisfied [62]. In this paper, we will let n = 2 for simplicity. 
Constant curvature solutions (for example: de Sitter solution: R = const = XiRq > 0) 



are the roots of the algebraic equation 



62| 



Rf{R)=2f{R). 



(6) 



Substituting the expression of f{R) given in Eq. ([5]) into the above equation, one can 
obtain 

' 2[{1 + xiy^+^ - 1 - {n + i)xiy 
In order to satisfy the stability conditions of the system, the following inequality must be 



satisfied 



62| 



(1 + > 1 + (n + 2)xf + (n + l)(2n + 



(8) 



Setting n = 2 and solving the above inequality, one gets Xi > y v 13 — 2, which leads to 
Xs > 0.94. We use = 0.95 in this paper, without loss of generality. 



III. THE GROWTH OF MATTER PERTURBATIONS 



As shown in Refs. 1171 . Ill8l |. the background evolution of a viable f{R) is very compli- 
cated. Here, we neglect all higher derivative and non-linear terms, and we then obtain the 

4 



equation governing the growth of matter perturbations on subhorizon scales as follows 99 1 



5 + 2H5-A-nGeffpra5 = 



(9) 



where Ggjf i s an effective Newton gravity constant and for an f{R) model, it can be 



expressed as 



103 | 



F 1 + 3 



k^F' ■ 
a'^F 



Defining the growth factor / = dln6/dlna, Eq. (Q becomes 



din a 2 \ 



(10) 



:ii) 



d\na J 2 Gjv 

Obvio usly, the growth factor is scale dependent, which leads to a dispersion of growth 
index 104| . Here we consider the wavenumber k within the range 



0.01 h Mpc-^ <k<0.2h Mpc-^ 



(12) 



which is relevant to the galaxy power spectrum 1161] . In scale smaller than 0.2 h Mpc~^, 
non-linear effects are obvious and for scale larger than 0.01 h Mpc~^ the current obser- 
vations are not so accurate. 



A. a constant 7 

In this subsection, we discuss the parameterized form / = ~ with a constant 
7. Usually, the theoretical value of 7 can be obtained by expanding the equation of 7 
around Qm — 1, which is a good approximation at the high redshift. In principle, we 
can also obtain the theoretical value of 7 by solving Eq. (11 II) numerically and using the 
value of Qmo given by current observations. Since the observational results on Q^no for 
Starobinsky's mode l is n ot obtained yet, we use Qmo = 0.278]'^o;g23 at the 68% confidence 



level given in Ref. |l05| with a model independent method. Solving Eq. ( ITTi) to obtain 
/(O) numerically and using the relation /(O) = fi^Q with Q^no taking the best fit value 
0.278, we find 70 — 0.42, which seems to be almost independent of the value of k. 

In order to discriminate different models with the growth factor, we must compare 
the theoretical value and the observational one of 7. The current observations give 12 



data poi nts o 



in Refs. 



112 



the growth factor |l06l - lll3| . Let us note that although the data given 



113j are measured without 'any' bias, other data points are obtained by 



assuming a flat ACDM model with Qm,o taking a specific value, for example, Qrn,o = 0.25 
or 0.30. So, caution must be exercised when using these data. With this caveat in mind. 



il 



114 



115| . Using these 12 

+0.17+0.47 



it may still be worthwhile to apply the data to fit models 

data, we find that, for a constant 70 and ilmo = 0.278, = 4.6 and 7 = 0.63!^o;;[4^o;|3 at 
the 1(7 and 2a confidence levels. It is easy to see that the Starobinsky's model is allowed 
only at the 2a confidence level. However, by comparing / and fij^, we can see that the 
error rate is larger than 10% as shown in Fig. ([1]), which means that the result obtained 
with a constant 7 may be biased. This bias arises from the fact that 7 is a function of 



redshift instead of a constant. More recently, the authors in Refs. 9ll-l95| discussed a 
linearized form of 7 with 7 = 70 + 71^, where 71 = 7g = ^(2; = 0). In the subsequent 
subsection, we will examine this varying form of 7 in detail. 



B. 7 = 70 + liz 

This linearized form of 7 has been studied in the wCDM, DGP and f{R) gravity, and 
it gives a very good approximation at the redshift region 2; < 0.5. In Ref. IqgI, we found 
that the constraints on 70 and 71 from three low redshift observational data cannot rule 
out the DGP model at la confidence level. Here, we want to see what happens for the 
Starobinsky's model, where we have 

71 = [ln^];,|o]" 

When Qm,o = 0.278, we obtain that 70 — 0.41. At the same time, we find that, for 
different k, the variation of 71 is small, for example, 71 varies from —0.20 to —0.24 when 
k is from k = 0.0.1 h Map-^ to 0.2 h Map-^. In Fig. 02]), we give the relative difference 
between the growth factor / and Q'^^'^^^ with Qmo = 0.278 and find that, at low redshifts, 
the error is below 2%, which means that this linearized form gives a better approximation. 

Now we discuss the constraints on 70 and 71 from the observations. Since this linearized 
form is valid in the low redshifts, only three low redshift data points can be used. Fig. (|3]) 
shows the results. From this figure, one can see again that only at the 2a confidence level 
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-^:^'o-3(7o-7.)(-l 



2Ho 



2 G 



N 



"mfl 



1 

2 



(13) 



is the Starobinsky's model consistent with the observations and it can be ruled out at the 



96|. 



la confidence level. This is in contrast with the DGP model 

However, the approximate form / ~ Q'^'^'^'^^ is only valid at the low redshifts. In order 
to use all the current observational data, we need to find a new approximate expression 
of /, which can give a good approximation in all redshift regions. 



C. a new approximation of / 



From the Fig. (2) in Ref. [9ll-l95|. which gives the evolution of /, one can see that / is 
larger than 1 in the region of 1 < z < 3.5. Since in a fiat universe Qm is always less than 
one, according to the usual approximation / ~ f^m^-* one cannot obtain /> 1 if 7(2;) is 
positive in the region 1 < z < 3.5. Thus the usual parameterized form of / is hard to 
give a good approximation. 

From the definition of Ge// and the fact that F is close to one for z > 1, one has at 
the redshift region z > 1 



Gn 1 + 3^ 



We assume an approximation of / by multiplying fij^ with a factor similar to the above 
expression, i.e., we assume 

where a is a constant. We find that when a is about equal to 0.85, for different wavenum- 
bers fc, the error rate, which is defined as (1 + a (i^z)'^^^ )^m/ f " 1; is about less than 
5%. The result is shown in Fig. Therefore, with this new parametrized form of /, 
one may use all the observational data. After numerical calculations, we obtain that 70 is 
about 0.57, which seems to be almost independent of the value of wavenumbers k, when 
a = 0.85. Using the 12 observational data points of the growth factor, we place the 
constraints on a and 70, which are shown in Fig. From this figure, we still find that 
the Starobinsky's model is allowed by the current observations only at the 2cr confidence 
level. 
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IV. CONCLUSIONS 



In this paper, we study the growth of matter perturbations in an f{R) model proposed 
by Starobinsky. Firstly, we discuss the case of a constant growth index. By comparing 
the theoretical value and the observational one, we find that the Starobinsky model is 
allowed by the current observations only at the 2a confidence level. However, in this 
case, the error rate between the growth factor / and is larger than 10%, so, the 
result obtained with a constant 7 may be biased. Then, a linear expansion of growth 
index, 7 = 70 + 712;, is studied, which is valid at the low redshift region z < 0.5 and 
gives a better approximation at these redshifts. With three low redshift observational 
data, we find again that the Starobinsky model is allowed only at the 2a confidence level. 
Finally, in order to use all the present data, we propose a new approximate form of /, and 
show that this new form gives a reasonable approximation both at low and high redshift 
regions. For different scales, the largest error is less than 5%. With this new proposed 
form of /, we still find that the Starobinsky model is consistent with the observations 
only at 2a confidence level. So, our results seem to suggest that although the Starobinsky 
f{R) model is excluded by the current growth factor data at la confidence level, it is still 
allowed at 2a level. 

It should be pointed out that, in our discussion of the growth of matter perturbations, 
the higher-derivative terms were di scard ed. Recently, it has been found, that with the 
covariant perturbation theory (see 119| for a recent review), which offers the simplest 
way to describe the evolution of the perturbations, these higher-derivatives terms can be 
kept in the analysis of matter growth. So, it remains an interesting topic to examine what 
happens when the effects of these terms are taken into account. 
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FIG. 1: The relative difference between the growth factor / and rimwith 0,^,0 = 0.278. 
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FIG. 4: The relative difference between the growth factor / and (l + a (i-\-zy^+s )^^ with Q, 
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FIG. 5: The la and 2a contours of 70 and a by fitting the Starobinsky's model with the current 
growth factor data. 
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